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We deduce the mass of the pseudoscalar glueball G from an r]-r]'-G mixing formalism based on the 
anomalous Ward identity for transition matrix elements. With the inputs from the recent KLOE 
experiment, we find a solution for the pseudoscalar glueball mass around (1.4 ± 0.1) GeV, which 
. is fairly insensitive to a range of inputs with or without Okubo-Zweig-Iizuka-rule violating effects. 

' This affirms that r;(1405), having a large production rate in the radiative J/^ decay and not seen in 

77 reactions, is indeed a leading candidate for the pseudoscalar glueball. Other relevant quantities 
C , including the anomaly and pseudoscalar density matrix elements are obtained. The decay widths 

for G 77, i~ are also predicted. 
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I ■ I. INTRODUCTION 

^: 

(— I ' The quest for pseudoscalar glueballs has continued for decades. An i?(1420) meson with a mass of 1426 MeV 

was first discovered at CERN in 1963 through pp interactions [l|. In 1980, Mark II observed that a J/ip meson 
decays via photon emission into a resonance at a mass around 1440 MeV This new state, named i(1440) 
by Mark II and Crystal Ball Collaborations [3j], was also once called G(1440) in Shortly after the Mark 

II experiment, i?(1420) and t(1440) were first proposed to be the pseudoscalar glueball candidates in Q and in 
respectively, while an opposite opinion that i?(1420) was an 1+ ss quark state was advocated in [8|. 
If) . As the experimental situation was sorted out, i?(1420) turned out to be an 1+ meson now known as /i(1420), 
CnJ ' and t(1440) was a pseudoscalar state now known as ?7(1405). For an excellent review of the E and l mesons, 
see p. 

77(1405) indeed behaves like a glueball in its productions and decays. The KKn and rjiTTT channels in 77 
00 ! collisions have been investigated ^u\. While ?7(1475) in KKtt was observed, 77(1405) in 777r7r was not. Since the 

■ glueball production is presumably suppressed in 77 collisions, the above observations suggest that the latter 
^ I state has a large glueball content [11|. J/t/; radiative decays through jgg have been considered as the ideal 

■ channels of searching for glueballs. The branching ratio B{J /ip 77^(1405)) of order 10"'^ is much larger than 
' the decays J/-0 — > 777(1295), 777(2225),... which are either not seen or are of order 10~*. The decay of a nearby 

77(1475) — ^ 77 has been observed [l3], but 77(1405) 77 has not. All these features support the proposal that 
77(1405) is a good pseudoscalar glueball candidate 0. There were also theoretical support based on the closed 
flux-tube model [l^l and the model that combines the octet, the singlet, and the glueball into a decuplet [13]. 
Besides 77(1405), other states with masses below 2 GeV have also been proposed as the candidates, such as 
77(1760) in [H and ^(1835) in [H. 

As for the scalar glueball, two of the authors (HYC and KFL) and Chua [l^ have considered a model for 
the glueball and qq mixing, which involves the neutral scalar mesons /o(1370), /o(1500) and /o(1710), based on 
two lattice results: (i) a much better SU(3) symmetry in the scalar sector than in the other meson sectors [l3| 
and (ii) an unmixed scalar glueball at about 1.7 GeV in the quenched approximation [l^. It was found that 
/o(1500) is a fairly pure octet, having very little mixing with the singlet and the glueball, while /o(1370) and 
/o(1710) are dominated by the glueball and the qq singlet, respectively, with about 10% mixing between them. 
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The observed enhancement of a;/o(1710) production over (/)/o(1710) in hadronic J/i/' decays and the copious 
/o(1710) production in radiative J /ip decays lend further support to the prominent glueball nature of /o(1710). 

Contrary to the above case, the pseudoscalar gluebaU interpretation for 77(1405) is, however, not favored by 
quenched lattice gauge calculations, which predicted the mass of the 0""*" state to be above 2 GeV in [l^ and 
around 2.6 GeV in [H, [l^l- It is not favored by the sum-rule analysis with predictions higher than 1.8 GeV 
[UHSl either. Readers are referred to 23] for a recent review on the results of the glueball masses. Note that the 
above lattice calculations were performed under the quenched approximation without the fermion determinants. 
It is believed that dynamical fermions may have a significant effect in the pseudoscalar channel, because they 
raise the singlet would-be- Goldstone boson mass from that of the pion to rj and 77'. It has been argued that the 
pseudoscalar glueball mass in full QCD is substantially lower than that in the quenched approximation |22l] . In 
view of the fact that the topological susceptibility is large {~ (191MeV)^) in the quenched approximation (2^ . 
and yet is zero for full QCD in the chiral limit, it is conceivable that full QCD has a large effect on the glueball 
as it does on rj and 77'. 

In this paper, we infer the pseudoscalar glueball mass mc from the rj-rj'-G mixing, where G denotes the 
physical pseudoscalar glueball. Implementing this mixing into the equations of motion for the anomalous Ward 
identity, that connects the vacuum to 77, 77' and G transition matrix elements of the divergence of axial- vector 
currents to those of pseudoscalar densities and the U(l) anomaly, ma is related to other phenomenological 
quantities such as the 77, rj' masses, the decay constants, and the mixing angles. Since the mixing angles have 
been measured recently from the IVtIV' decays by KLOE (25j . mo can be solved. Our numerical study 
gives a fairly robust result mc « 1.4 GeV, which is insensitive to a range of inputs. We also obtain the matrix 
elements for the pseudoscalar densities and axial U(l) anomaly associated with the 77, r]' , and G states. The 
values of the pseudoscalar density matrix elements for the 77, 77' mesons are close to those obtained in the 
Feldmann-KroU-Stech (FKS) scheme [2^, which does not include the mixing with the pseudoscalar glueball. 
The results of the anomaly matrix elements for the above states are quite consistent with those estimated 
from the topological susceptibility [13, [H, IIH and the lattice evaluation [li| , indicating that the J/ip ^ jrj' 
branching ratio could be comparable to that of Jf-ip 7G. We then study the pseudoscalar glueball decays 
into two photons and two leptons G — s- 77, £^£^ . The comparison of our analysis with the properties of known 
mesons suggests that the 7/(1405) meson is a strong pseudoscalar glueball candidate. 

In sec. II we set up the formalism for the rj-rj'-G mixing, assuming that the glueball only mixes with the 
flavor-singlet 771, but not with the flavor-octet 778 . Our parametrization for the mixing matrix contains only 
two angles and differs from that in [255, where it is assumed that rj does not mix with the glueball state. The 
solution for the pseudoscalar gluaball mass 77^0 is derived in Sec. Ill with the phenomenological inputs from 
KLOE [2^. The solutions with the inputs from [s^l and from [2^ as a limit of vanishing mixture with the 
glueball state are also presented for comparison. It will be shown in Sec. IV that the result for rric is stable 
against the variations of phenomenological inputs and of corrections violating the Okubo-Zweig-Iizuka (OZI) 
rule [3l|. The G 77, decay widths are also estimated. Section V is the conclusion. 



II. Tj-ri'-G MIXING 

We extend the FKS formalism [2^ for the 77-77' mixing to include the pseudoscalar glueball G. In the FKS 
scheme, the conventional singlet-octet basis and the quark-flavor basis have been proposed. For the latter, the 
qq = {uu -\- dd)/\/2 and ss flavor states, labeled by the 77, and 77s mesons, respectively, are defined. In the 
extension to the rj-rj'-G mixing, the physical states rj, rj' and G are related to the octet, singlet, and unmixed 
glueball states 7/8, 771 and 5, respectively, through the combination of rotations 



with the matrices 



umuiic^a) \m) , (1) 





1 

U3{9)= I sine cos6' | , f7i(0G) = | cos0g sin0G | ■ (2) 

— sin cos (pG 
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The matrix Ui {U3) represents a rotation around the axis along the rjs meson (unmixed glueball g). Equation ([T]) 
is based on the assumption that rjs does not mix with the giucbaU, under which two mixing angles and (j)G 
are sufficient. 

The octet and singlet states are related to the flavor states via 

/ M \ 

= C/3(0.) \ris) , (3) 
V \9) J 

where 9i is the ideal mixing angle with cos6i — \/T/3 and sin^i = ^2/3, i.e., 9i — 54.7°. The flavor states are 
then transformed into the physical states through the mixing matrix 

cos (J) + sin 9 sin 0i — sin (j) + sin 9 cos 0^ Ac — sin 9 sin (pa 
smcj) — cos9 siii9iAG cos — cos 61 cos 6',; Ag cos^sin^c | , (4) 
— sin 9i sin — cos 9i sin cos 

with the angle (p = 9 + 9i and the abbreviation Aq = 1 — cos (po ■ U has been written in the form, which 
approaches the FKS mixing matrix ^| in the (?!)g — > limit. That is, the angle (f> plays the same role as the 
mixing angle in the FKS scheme. 

Our formalism assumes isospin symmetry, i.e. no mixing with tt*', and neglects other possible admixtures 
from cc states and radial excitations. The widely studied decay constants fq and fs are defined by (26j 

{Q\sri5s\7ls[P)) - -ifsP^ , (5) 

for the light quark q = u 01 d. The ijq (77^) meson decay constant /| (/|) through the s {q) quark current j32i ]. 
and the unmixed glueball decay constants /J^'" through the q and s quark currents, can be defined in a similar 
way: 

{0\qri5q\Vs{P),9{P)) = -^fljP'' ' 

(0|s-7^5s|r/,(P), g{P)) = -^flq P^ . (6) 
The decay constants associated with the 77 meson, rj' meson, and the physical glueball defined in 

mi^l,qHP).v'{P)MP)) = -^fl^'^aP" > 

(0|s-7^75s|r?(P), 77'(P), G(P)> = -i.f^,,',G P" , (7) 
are related to those associated with the rjq, ijs, and g states via the same mixing matrix 

fn fri \ j fq fq \ 

4 fv' ^U{cf>,<l>G) f! fs . (8) 
fh fhl \fl ft) 

Sandwiching the equations of motion for the anomalous Ward identity 

dMl'^l^q) = 2imq q-f^q + G^. G^'' , 

47r 

a^(s7^75s) = 2zm, + ^ G^, G^^^ , (9) 

47r 

between vacuum and \rf), jry') and |G), where G^^ is the field-strength tensor and G^'^ the dual field-strength 
tensor, and following the procedure in [33|, we derive 

= C/t(0, <j>G)M^UicP, ^g)J . (10) 
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In the above expression the matrices are written as 

/ ml, + (^/2//,)(0|a,GG/(47^)|,?,) m% + (l//,)(0|a,GG/(4^)|77,) 
Ml, = I m^^ + (y2//,)(0|a,GG/(47r)|,7,) + (l//s)(0|a.GG/(4^)|77,) 

{V2/f,masGG/{47r)\g) m^, + (l//.)(0|a.GG/(4^)|g) 
1 





= \ m^, I , J = ( f!/f, 1 I , (11) 

/g / fq fg/ fs 



with the abbreviation 

a/2 

mlq,qs,qg = {0\mumj5U + mddi'^5d\r]q , Tjs , g) , 
Jq 
2 

= -^(0|m.ssi75s|779,?7s,ff) • (12) 

Is 

In the hmit of the large color number N^^ the scaling for the decay constants, the pseudoscalar densities, and 
the anomaly matrix elements is [3^ 



fq,s ^ 






-0(1), 


fg - n - 0(l/y7Ve) , 


TOG ^ 


^O(l), 


4>G ^ 




1 


™gq ' 


-o(i), 




-0(1), 




mlg^ 


^ ml, ^ 0(l/\ 






m% ^ ml, ^ 0{l/N,) , 



(0|a,GG/(4^)|5) ^ 0(1), (0|a,GG/(47r)|77,) ^ (0|a,GG/(47r)|77,) ^ 0(l/yA^) . (13) 

The pseudoscalar meson and glueball masses scale as 0(1) in large Nc- However, it has been pointed out 
' ' ' ' 2 



[271, |28|, |29| that the sub-leading 0{1/Nc) term in the 771 mass squared to^^ ~ 0(1) + 0{1/Nc) is numerically 
large due to the U{1) anomaly, and is related to the topological susceptibility x in the quenched QCD without 
fermions. In the chiral hmit, the relation m^^, = ANpx/fl = 2VA^F(0|asGG/(47r)|?7') //^r with Np being the 

number of flavors gives (0|asGG/(47r)|?7') w 0.035 GeV^ for x = (191 MeV)" H. Although being 0(1//^) 
in large Nc, this anomaly matrix element is numerically much larger than the 0{\/ Nc) quantities ml,f, « 

0.0026 GeV'^ for ml, « toJ and comparable to ml^fs « 0.087 GeV"^ for mlg ~ 2m\ — toJ. In view of this, we 
shall keep all the anomaly matrix elements for r]q,r]s and g in the following analysis. On the other hand, we 
expect the decay constants and the pseudoscalar density matrix elements to have the normal ordering in terms 
of Nc- That is, we expect fq^s > fg'^ > f,, fsi m?,, > ml, > to^^, and m^^ > ml, > ml,. The above ordering is 
consistent with the OZI rule in that double quark annihilation, which is present in /|, fj, ml^ and ml, but not 
in others, is OZI-rule violating and suppressed. We note that the two sides of each of the equations in Eq. (jTU]) 
have the same Nc scaling, implying consistency of our formalism in terms of Nc- 
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III. PSEUDOSCALAR GLUEBALL MASS 



The exphcit expansion of Eq. pO| leads to 

U^Un + U12R' + [/i3r']m' + UUU21 + U22R' + C/23r']m^' + uUUai + U32R' + Ussr'jml 

= m% + (V2//,)(0|a,GG/(4^)|77,) , (14) 
UliUiiR + U12 + Ui3r]ml + UUU21R + U22 + ^723^]™^, + U^p^iR + U32 + f/33?-]m^ 

= m% + (l//,)(0|a,GG/(4^)|77,) , (15) 
= ml + (V2//,)(0|a,GG/(47r)h,) , (16) 
[/2\[C/iii? + U12 + Uur]ml + UUU21R + U22 + U2zr]ml, + U^U^^iR + C/32 + Ui3r]ml 

= mL + (l//.)(0|a.GG/(47r)|r7,) , (17) 
UliWii + U12R' + C/i3rV' + Ul^[U2i + U22R' + C/23'^V'' + Ul^[Uzi + U^2R' + U^zAml 

= m^^ + (V2//,)(0|a,GG/(4^)|g) , (18) 
Ul^piiR + U12 + Ui^r]ml + [/j2[(72ii? + U22 + t/23r]m2, + U^U^iR + [732 + !733r]m2j 

(l//,)(0|a,GG/(4^)|g) , (19) 



TO 



,2 



where the parameters r = fg/fs, r' = f^/fq, R = fq/fs, and R' = f^/fq are introduced, and Uij denotes 
the matrix element of U. In developing our mixing formalism, the flavor-independent couplings between the 
glueball g and the pseudoscalar mm, dd and ss states are assumed, so that g only mixes with the flavor-singlet 
rji . Accordingly, we postulate = V^fg and /I — fq, and thus the relations 

r' = V2j-r R' = ^i?, (20) 
Jq Jq 

which will be adopted in the numerical study in Sec. IV. 

We first explore the implication of the rj-rj'-G mixing on the glueball mass toq. To simplify the matter, the 
ratios r, r', R and R' are neglected, which are 0(l/\/^) and 0{1/Nc), respectively, in large Nc as shown in 
Eq. We also neglect m^^ and to^^ relative to the numerically large anomaly term {Q\asGG / {ATT)\g) / fqg as 

an approximation. It should be safe to drop m^g, since it is, like the small to^^ « m^, proportional to the light 
u/d quark mass. On the other hand, it is not clear if it is safe to drop to^^. Although it is 0(l/\/7Vc) compared 
to m^s, but the latter, being proportional to the strange quark mass, is larger than {Q\asGG / {^'K)\r]') / f ^ in the 
chiral limit as discussed in Sec. |lTl This subject will be investigated in a more detailed numerical analysis later 
in Sec. IIVI Having made the above assumptions, we take the ratio of Eqs. and and obtain 

ce{scj) - ces0,l^G)ml, - se{c(j) + sese, AG)^ml - se^ccj)Gm% _ V2/^ 
ce{c(j) - cec9^/:\G)ml, + s9{s(f) - sece.Aayml - c0iC(l)Gm% " fq ' 

where c(j) (s0) is the shorthand notation for cos (fi (sin (p) and similarly for others. 

Note that the above simple formula still holds, even after keeping the r'- and r-dependent terms, as long as 
they obey Eq. (PO)) . In other words, the r'-dependent terms in Eq. ^TE\\ and the r-dependent terms in Eq. 
can be absorbed into the right-hand sides of these equations and are therefore canceled after taking the ratio 
of Eqs. (|18p and (|19p. The factor sin^c in the numerator and the denominator of Eq. (^J) has been canceled 
out, so that the 0g dependence appears at order of Aq « 0^ for small (pG- As such, we find that the solution 
for TUG is stable against the most uncertain input 0g, as long as the 77,77' mesons do not mix with the glueball 
g intensively. The solution depends on the ratio fs/ fq, which is, through Eqs. and (|16p . related to 

V^fs ^ (0|a.GG/(4^)|7?,) 
fq (0|a,GG/(47r)h,) ' ^ 

if the pseudoscalar density matrix elements m^^ and to^^ are neglected. It implies that the SU(3) symmetry 
breaking in the axial anomaly matrix element echoes the symmetry breaking in the decay constants, and plays 
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a sensitive role in the determination of the pseudoscalar gluebaU mass. For a given (pc, ma increases with 
decreasing fs/ fq- By the same token, when the anomaly matrix element for the rjs becomes larger relative to 
that of the r]q meson, the mass of the pseudoscalar glueball gets higher. 

Before solving for mc from Eq. pT|) . we explain the strategies for data fitting adopted in [IE] and [13], which 
led to different extractions of the mixing angle 0g. In Ref. [1^, the decay constants fq = (1 ± 0.01)/^ and 
fs = (1.4 ± 0.014)/^ [s^l, and the parameters associated with meson wave function overlaps [3^ were fixed as 
inputs. The angles (j> = (39.7 ± 0.7)° and (j)G = (22 ± 3)° were then determined from the relevant data. A tiny 
error (1%) was assigned to fq and fs, which is one of the reasons why a high precision was reached for the 
determination of 0g- In the data of P ^ and V were first considered, which do not depend on fq 

and fs, and the fit gave the outcomes (j> = (41.4±1.3)° and = (12±13)°. Without precise inputs of fq and fs, 
and with the hadronic parameters for meson wave function overlaps being free, it is not unexpected to get a wide 
range for (f>G- The value of being consistent with zero [1^, means that the data could be accommodated 
by the hadronic uncertainty alone. The extracted </> and 0g were then used as inputs to determine fq and fs 
from the r/, ?/ — > 77 data. Since 0g has a wide range, the results fq = (1.05 ± 0.03)/,r and fs = (1.57 ± 0.28)/7r 
also have larger errors. The correlation between (f>G and fs (a smaller 0g corresponding to a larger fs) is a 
consequence of the constraint from these data. We also note that a larger mixing angle (j)G — (33 ± 13)° has 
been extracted from the J/ip — VP data recently . In summary, both sets of parameters in [s^] can fit 
the data, and are consistent with each other within their uncertainties. It is seen that fq, fs, and (f> are more 
or less certain, but 0g varies in a wider range. Fortunately, the solution for the pseudoscalar glueball mass ttiq 
is not sensitive to 0g as discussed above and will be explored further in the remainder of this paper. 

As stated before, KLOE postulated that the glueball does not mix with 77 |^]. We shall point out that this 
postulate does not yield a solution for ma in our formalism. The KLOE parametrization for the rj-rj'-G mixing 
matrix is written as 

(cos 4> — sin (j) 

sin (f) cos 4>G COS0COS0G sin^G | . (23) 
— sin (j) sin (f)G — cos (j) sin cf)G cos (f)G 

Repeating the above procedure, Eq. (|2ip is modified to 

s(j)C(l}GS(j)G{ml, ~ m%) s(f> ^2/, 



c(t)c4)GS(j)G{rn^^, ~ my) c(j) fq 



(24) 



For the KLOE parameter set fq = /^r, fs — ^■'^f-r: and = 39.7° 25], Eq. ([24]) does not hold, and there is no 
solution for tug as a result. 

Since we have changed the mixing matrix from KLOE's in Eq. to Eq. ([U, we need to refit </) and (t>G in 
principle. However, comparing Eqs. and ([3]), it is easy to find that their 2x2 sub- matrices in the left-upper 
hand corner have almost equal elements for « 40° and (/)g ~ 22°: 

/ 0.751 -0.654 0.097 \ / 0.766 -0.643 \ 

U = 0.585 0.725 0.362 , C/rloe = 0.596 0.710 0.375 . (25) 
\ -0.306 -0.216 0.927 / \ -0.241 -0.287 0.927 / 

These four elements, which are responsible for the quark mixing, are the only ones involved in the data fitting of 
4> jr],jri', rj' — > 7/5,70;, and 77,77' — > 77 mentioned above. Therefore, it is expected that the refit of the data 
using our parametrization will give the mixing angles close to KLOE's. That is, the KLOE parameter set can 
be employed directly in our numerical analysis within uncertainty. The other off-diagonal elements in Eq. (|25p , 
describing the mixing among the 7/, ?/ mesons and the glueball, do have different values. It is thus understood 
why the two parametrizations have similar mixing angles, but the ratios in Eqs. (PT|) and ([M)) exhibit different 
behaviors as far as ttt-g is concerned. 

It is also interesting to consider the parameter set from |26|] with fq — (1.07 ± 0.02)/^, fs — (1.34 ± 0.06) fn, 
(j) = (39.3 ± 1.0)°, and i^g = (no mixing with the pseudoscalar glueball). Note that the lower fs in [20| 
arises from combined experimental and phenomenological constraints. If only the experimental constraints were 
used, mainly those of the 77, 77' 77 data, its central value would increase and the range is enlarged, gi ving 
fs = (1.42 ± 0.16)/,r close to that extracted in [25]. Using the central values of fs/fq and (pG from |25l. l2a. |30|] 
as inputs, we derive the pseudoscalar glueball mass from Eq. ([2T|) [see also Eqs. (|27 |) - (f29| below] 

niG = 1.41, 1.56, 1.30 GeV , (26) 
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respectively. The above investigation leads to ma = (1.4±0.1) GeV with the currently determined phenomeno- 
logical parameters. The proximity of the predicted mc to the mass of 77(1405) and other properties of 77(1405) 
make it a very strong candidate for the pseudoscalar glueball. We shall come back to visit the robustness of 
our prediction in the next section, when higher order effects in 1/Nc are included . 

One may question whether other pseudoscalar mesons with masses around 1.4 GeV, such as 7^(1295) and 
7/(1475), should be included into our mixing formalism. We note that 7/(1295) and 7/(1475) have been assigned 
as the 2^5*0 states, namely, the radial excitations of qq and ss, respectively [sl, [3§|. As stated in the previous 
section, these radial excitations do not mix with the rj and rj' mesons by definition, since they are diagonalized 
under the same Hamiltonian. As for the mixing of the radial excitations with the glueball, we speculate that it 
is negligible for the following reason. 77(1295) is practically degenerate with the radially excited pion 7r(1300), 
and 7/(1475) is about 200 MeV above 7/(1295), a situation similar to the ideal mixing in the vector meson 
sector with 0(1020) being 200 MeV above w(780). This suggests that 7/(1295) and 77(1475) are much like the 
radially excited isovector pseudoscalar qq and pseudoscalar ss without annihilation. The difference between the 
isoscalar 77, 7/' mesons and the isovector pion is that the former have disconnected insertions (annihilation) with 
the coupling going through the contact term in the topological susceptibility which pushes their masses up. By 
virtue of the fact that 77(1295) is degenerate with 7r(1300) and 7/(1475) is ~ 200 MeV above, they do not seem to 
acquire such an enhancement for their masses. Therefore, we venture to suggest that the annihilation process 
is not important for these two mesons, and their mixing with the glueball is weak. 

A pseudoscalar glueball mass about 1.4 GeV was also determined from the framework of the rj-rj'-G mixing 
in [40], but with a strategy quite different from ours: The mixing is assumed to occur through a perturbative 
potential, so that the mixing angles are parametrized in terms of the transition strength among the states 7/g, 
77s and g and their masses 771^^ , m,^^ and rrig 40] . These parameters were then fixed from data fitting. Hence, it 
is the unmixed glueball mass nig, instead of the physical glueball mass trq, that was derived in [40j. Moreover, 
the result of ^40/] is a consequence of data fitting, while ours comes from the solution to Eq. (|2ip . If the quark 
flavor states do not mix strongly with the glueball, 777^ is expected to be close to that of 777^. Following this 
reasoning, three possible glueball candidates, 7/(1405), 7/(1475), and X(1835) with masses around 1.4 GeV, 
have been examined in [40j . and the latter two were found to be experimentally disfavored. 

IV. NUMERICAL ANALYSIS 

We now proceed to solve Eqs. p^ - P^ based on the large Nc hierarchy in Eq. As discussed in Sec. HIl 

all the anomaly matrix elements will be kept. Even though they are small parametrically (0(1) and 0{l/ \/Nc)), 
they are large numerically. As a follow up of the last section, we first neglect the decay constants /|, /|, and 
/|'*, which are 0{l/y/Nc) and 0{1/Nc) lower than fq,s, respectively. We also neglect the pseudoscalar density 
matrix elements rn'^g, m^g, nigg, and 771^^, which are similarly suppressed as compared to 771^^ and 777^,,-^. Under 
this approximation, our formalism involves six unknowns: three mass related terms ttt-g, w^^, and 777^^, and 
three anomaly matrix elements (0]asGG/(47r)]7/g), (0]asGG/(47r)]7/s), and {0\asGG / {4:Tr)\g) , provided that the 
phenomenological quantities 777^, 777^,, fq, fs, 4> and 0g are given as inputs. There are six equations from 
Eq. (jlOp . so the six unknowns can be solved in principle. We note in passing that the four unknowns 777^^, 777,^^, 
(0]asGG/(47r)]7/g), and (0]Q;sGG/(47r)]7/s) were solved for five given inputs 777^^, 777^,, fq, fs, and <j) in the rj-rj' 
mixing case '41] . 

Using the central values of the parameter sets from [2^, [sO], and [2^ for fq,fs,4> aiid 4>G as inputs, we obtain 
the following solutions 

= -0.073 GeV^ , (OlasGG/(47r)l77q) = 0.069 GeV^ , 
7772^ = 0.52 GeV^ , (OlasGG/(47r)l77s) = 0.035 GeV^ , 

777G = 1.41 GeV , (0lasGG/(47r)lg) = -0.033 GeV^ , (27) 



^ The off-diagonal mass terms m^^ and m^^ have been absorbed into the matrix elements (0|asGG/(47r)|7jq) and (0|«s'-'G/(47r) |r;s ), 
respectively in |26|| . 
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TABLE I: Solutions for various inputs of (0|aaGG/(47r)|77,) = 0.065 GeV^ (the first row), 0.050 GeV^ (the second row), 
and 0.035 GeV'^ (the third row) with niqq = m^, — 2m^ — m^, rrisg — m^g = m^s — rrigq — and (j> — 42.4°. The 
upper (lower) table is for (j)G = 22° {(j)G = 12°). 



r 


R fs 


mG(GeV) (0| 


f^GGM (GeV^) 


(OmCGlg) (GeV^) 


0.004 


-0.002 1.25/^ 


1.50 


0.037 


-0.038 


0.22 


-0.002 1.25/^ 


1.50 


0.028 


0.036 


0.44 


-0.002 1.25/^ 


1.50 


0.020 


0.111 


-0.26 


-0.003 1.28/^ 


1.44 


0.036 


-0.108 


0.16 


-0.003 1.28/^ 


1.44 


0.028 


0.035 


0.58 


-0.003 1.28/^ 


1.44 


0.019 


0.178 



m^^ = -0.084 GeV^ , (0|a^GG'/(47r)|?7,) = 0.067 GeV^ , 

to2^ = 0.50 GeV^ , (0|a^GG/(47r)|?7,) = 0.032 GeV^ , 

mc = 1.30 GeV , (0|a^GG/(47r)|5) = -0.015 GeV^ , (28) 

and 

m^^ 0.012 GeV^ , (0|asGG/(47r)|r/^) = 0.054 GeV^ , 

ml^ = 0.50 GeV^ , (0|asGG/(47r)|77,) = 0.030 GeV^ , 

TOG = 1-56 GeV , (0|asGG/(47r)|g) = GeV^ . (29) 

The above solutions give us an idea of the range of uncertainties in our predictions. It is observed that the 
solutions for the anomaly matrix elements associated with the rjg and rjs mesons change little in Eqs. (f27|) - (f29|) . 
However, (0|Q;<;GG/(47r)|g) for the pseudoscalar glueball varies with the inputs oi (pc and fs/fq strongly as can 
be seen from Eqs. (fT8|) and (I19p . The solutions of to^^ in Eqs. ([27]) and ((28|) deviate from the naive expectation 
TO^q = toJ [2^ , while that in Eq. ([29]) is in better agreement with to^ due to a smaller fs . The solutions of m^^ , 
on the other hand, are stable with respect to the various inputs, and are close to the expected leading Nc result 
TOjg = 2to^ — TO^ . The values for tog have been shown in Eq. (|26p already. We should comment that to^^ is 
small because of the cancellation of two large terms as pointed out in [41J . It flips sign easily, depending on the 
inputs of fs/fq and OZI-rule violating effects, which have been considered before in the two-angle formalism for 
the rj-r]' mixing [13, ll^]. Our opinion is that introducing the OZI-rule suppressed decay constants J^, /J* is 
more transparent than employing the multiple- angle formalism. It has been observed that the tiny corrections 
from /| and can turn a negative to^^ into a positive value easily due to the smallness of m^^ |32l | . 

In the following, we investigate the higher 0{1/Nc) effects from the decay constants /I'", /|, and /I, i.e., 
from r and R on our solutions. The pseudoscalar density matrix elements to^^, mlg, nigs, and to^^ are still 
ignored. We take m^g = ml, mis = - ml and (0|asGG/(47r)|77,) = 0.065, 0.050, and 0.035 GeV^ (the 

typical values from Eqs. (f27|)-([29l)). as the inputs in order to solve for the unknowns r, R and fs- The relation 
m^s = 2to^ — ml seems to hold well for the earlier solutions in Eqs. (PTjl - iP^ . Thus, it is reasonable to fix it to 
its leading Nc value. Taking fg — f^, (jy — 42.4° and 0g = 22° and 12°, the corresponding solutions are listed 
in TablelH The results of i?, /s, and tog are independent of the inputs of (0|asGG/(47r)|77g), reaffirming that 
Tog is independent of r as seen from Eq. (|21[) . The magnitude of R is smaller than that of r, which in turn is 
smaller than unity. This finding is in agreement with the large counting rule. The decay constant fs turns 
out to be lower than those in [IS, [s^ , following from the observation that a smaller fs leads to a positive to^^ 
(4ll |. The values of tog and (0|QfsGG/(47r)|77s) are consistent with the range derived in Eqs. ((27ll - ((29|) . implying 
that these higher 0(1/A*'c) effects are small. The parameters r and (0|asGG/(47r)|(7) are found to be sensitive 
to the inputs, and both of them increase with decreasing (0|asGG/(47r)|?7g). 

Finally, we explore the impact of m^g on our solutions. To do so, we add fs as an input so that to^^ can be 
introduced as an unknown, to^^ is not considered, because its effect should be very minor as explained before. 
The results for the various inputs of fs = (1.24-1.30)/^, ipG = 22° and 12°, and {Q\asGG / {AiT)\rig) = 0.050 
(0.035) GeV'^ are listed in Table |TT] (jlTT}. In the large Nc analysis for the resolution of the U{1) anomaly 
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TABLE II: Same as Table |I] except that {0|asGG'/(47r)|r7q) = 0.050 GeV^ and fs is fixed to trade for m^g as a free 
parameter. 



/. r 


R 




mG(GeV) (0| 




1.24/. 0.22 


-0.001 


-0.009 


1.60 


0.028 


0.036 


1.26/. 0.22 


-0.003 


0.004 


1.47 


0.028 


0.036 


1.28/. 0.23 


-0.005 


0.016 


1.34 


0.028 


0.038 


1.30/. 0.24 


-0.007 


0.029 


1.21 


0.028 


0.040 


1.24/. 0.12 


0.001 


-0.054 


2.15 


0.027 


0.030 


1.26/. 0.13 


-0.001 


-0.029 


1.84 


0.027 


0.031 


1.28/. 0.15 


-0.003 


-0.005 


1.52 


0.027 


0.034 


1.30/. 0.24 


-0.005 


0.018 


1.16 


0.028 


0.045 



TABLE III: Same as Table |IT] except (0|a,GG/(47r)|r7,) = 0.035 GeV^ 



/. r 


R 


m% (GeV^) mc (GeV) { 


0|ffGG|r,.)(GeV^) (01 


^GG\g) (GeV^) 


1.24/. 0.40 


-0.001 


-0.009 


1.60 


0.019 




0.105 


1.26/. 0.45 


-0.003 


0.004 


1.47 


0.020 




0.113 


1.28/. 0.54 


-0.005 


0.016 


1.34 


0.021 




0.126 


1.30/. 0.69 


-0.007 


0.029 


1.21 


0.022 




0.148 


1.24/. 0.27 


0.001 


-0.054 


2.15 


0.018 




0.136 


1.26/. 0.34 


-0.001 


-0.029 


1.84 


0.018 




0.146 


1.28/. 0.51 


-0.003 


-0.005 


1.52 


0.019 




0.168 


1.30/. 1.18 


-0.005 


0.018 


1.16 


0.023 




0.262 



[2 71 128| . |29| , it is the combined contribution from a contact term and the glueball that cancels the 77' contribution 
to give a zero topological susceptibility in full QCD in the chiral limit. This combined contribution is just the 
topological susceptibility Xqucnch in the quenched QCD, which leads to the Witten-Veneziano mass formula 
m^, = 47Vi?Xquonch//T- Xqucnch IS Calculated to be « 0.00133 GeV* [2^, and the quenched glueball contributes 
about —11% to Xqucnch [l^ , which makes the contact term to be « 0.00148 GeV*. It is observed that the anomaly 
matrix element (0|asG'G'/(47r)|g) is 0.105 GeV^ or lar ger in Table [TTTl This anomaly matrix element contributes 
(0|a,.GG/(47r)|5)V(-4m2) w -0.00141 GeV* to the topological susceptibility for (0|a,.GG/(47r)|5) = 0.105 
GeV'^ and nig = 1.4 GeV. Namely, the glueball contribution is as large as but destructive to the contact term. 
As the glueball contribution and the contact term already cancel each other to a large extent, there is no room 
left for the contact term to cancel the sizable 77 and 77' contributions in order to end in a very small topological 
susceptibility in fuU QCD, which has a value x(fullQCD) = -(V5'0)/(l/m„ + l/nid + l/rris) - 4 x lO'^ GeV* 
[43 |. It implies that the anomaly matrix element (0|asGG/(47r)|(7) = 0.105 GeV'^ is probably too large. 

Based on the above reasoning, we believe that (0|Q;iGG/(47r)|(7) > 0.105 GeV'^ is not likely to be a viable 
solution. This criterion would exclude all the results in Table Hill with (0|Q;sGG/(47r)|77g) — 0.035 GeV'^ as an 
input. For Table HIl it is seen that m'j.g and niG do depend on fg sensitively. In some cases, we have toq as 
large as 1.84 GeV and 2.15 GeV, for which rn^g are negative and large. We cannot discard these solutions of 
rrigg a priori, but they are not favored due to their negative values. This issue can be sorted out, when lattice 
calculations of m^g with dynamical fermions are available. As fs > 1.30 /jr, mo becomes smaller than 1.2 
GeV, where there are no pseudoscalar glueball candidates. Therefore, if excluding the solutions with large and 
negative rri^g, the range (1.4 ± 0.1) GeV of the pseudoscalar glueball mass obtained in Sec. IIIII will be more or 
less respected. 

Having studied the higher 0(l/iVc) effects and confirmed that they are small, modulo the uncertainty regard- 
ing m^g, we shall simply use the typical results in Eq. [Eq. ([^5]) ] to obtain the anomaly matrix elements 
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for the physical states 77, 77' and G: 



3|a3GG/(47r)|77) 
'\a,GG/{4TT)\7j') 
^\a,GG/{4TT)\G) 



0.026(0.028) GeV^ , 
0.054(0.057) GeV^ , 
-0.059(-0.041) GeV^ 



(30) 



It is found that the value of (0|a^GG/(47r)|?7) is close to y/3/2f^ml/3 w 0.021 GeV^ obtained in and 
(0|asGG/(47r)|ry') is within a factor of two from its chiral limit estimated from the topological susceptibility, i.e., 
(0|asGG/(47r)|77') = 2y/NFx/f^ = 0.035 GeV^ Equation ^ also reveals that (0|asGG/(47r)|G) is almost the 
same as that from the quenched lattice QGD calculation, which gives |(0 1 GG/ (47r) | G) | — (0.06 ± 0.01) GeV'^ 
The fact that (0|asGG/(47r)|ry') is comparable to (0|asGG/(47r)|G), which defies the large Nc scaling in 
Eq. (|13p . implies that the rj' meson production in the J/ip radiative decay may have a branching ratio as large 
as that for the pseudoscalar glueball production. 

Given the mixing angles, we can predict the widths of the G 77, £^£^ decays, assuming that they take 
place through the quark content [III ■ The ratio of the G ^ 77 width over the 7r° 77 one is expressed as 



r(G ^ 77) 

r(7r'' —>■ 77) 



5-y^ sin^ 

Jq 



i sm (pa 



Js 



I sin (pG 



(31) 



We have confirmed that both the parameter sets in [25| and [30| give the 77, 77' — > 77 widths in agreement 
with the data r(77 77) « 0.51 keV and Ffr?' — * 77) ~ 4.28 keV [sll, by considering the similar ratios 
for the 77,77' mesons. The parameter set in [111 ([s^l) leads to a ratio 387 (83.3) in Eq. ((3T|) . i.e., the decay 
width r(G ^ 77) = 3(0.6) keV for r(7rO ^ 77) = 7.7 eV If 77(1405) is a pseudoscalar glueball, we 

predict the branching ratio S(77(1405) — > 77) = 6(1) x 10~^, i.e., an order of 10~^ for the total decay width 
r(77(1405)) = 51 MeV [11] . The above result can be confronted with future experimental data. The "stickiness" 
S has been proposed to be a useful quantity for identifying a glueball rich state [i^, which is defined as the 
ratio of T{J/ip 7G) to r(G 77) with the phase space factors taken out. Combining our predictions for 
the pseudoscalar glueball production and decay, we obtain S — 18-80 for G, which is much larger than S = 1 
as defined for the 77 meson. 

For the G £^£^ decays, we calibrate their widths using the available 7r° e+e^ and 77 IJ-'^I^^ data: 



r(G 



r(7rO 



r(G^MV-) 



) 



I f iriG 



9 \ 771^0 
3 



mG_ 
n 



5^ sin di sin <j)G 

Jq 



cos 6i sin cpG 

Js 



5 ^ sin 9i sin (pG 
m„ J \ Jq 



\f2^ cos di sin 



5-^ (cos (/) + sin B sin ft; Aq) — \Pi^ (sin + sin d cos Oi^^g) 
. Jq JS 



-2 



(32) 



For r(7r° e+e") = 4.8 x 10"'^ eV |38l, we obtain r(G e+e") = 1.9(0.4) x lO""^ eV using the parameter 
set from [H] ((H). For r(?7 ^ M+A*") = 7.5 x lO^^ eV [H, we have r(G ^ /i+A*") = 4.0(1.0) x 10"^ cV. 
If 77(1405) is a pseudoscalar glueball, the above predictions correspond to the branching ratios S(?7(1405) — > 
e+e") = 4(0.8) x 10"^^ and ;B(7;(1405) ^J-~^^^~) = 8(2) x 10"^*^, which would be quite a challenge to observe 
experimentally. 



V. CONCLUSION 



In this paper, we have formulated the rj-rj'-G mixing scheme via the vacuum to meson transition matrix 
elements for the anomalous Ward identity. The extension to include the glueball mixing with the flavor-singlet 
771 is a generalization of the FKS scheme for the 77-77' mixing [2^ . Therefore, only one extra angle (j)G for the 
mixing of the glueball state g and 771 is introduced in addition to the angle (j) in the FKS scheme. We have 
explained the different parameter extractions from the same set of 7y, 77' meson data in [25] and [30"], which give 
an idea of the uncertainties contained in the inputs. The obtained pseudoscalar glueball mass 771^ around 1.4 
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GeV is much lower than the results from quenched lattice QCD (> 2.0 GeV). It has been examined that our 
solution for niQ depends weakly on the ratio of the decay constants fs/fq in the favored phenomenological range 
and is stable against the variation of (pa and the higher 0{1/Nc) corrections. 

There may not exist a unique feature which tells a glueball apart from a quark-antiquark state. We need to 
combine information from J/Tp radiative decays, hadronic decays, as well as 77 and leptonic decays as advocated 
in [43|. The comparison of our solutions with the available data suggests that 77(1405), which is copiously 
produced in the J/ tp radiative decay but has not been seen in the 77 reaction, is a strong pseudoscalar glueball 
candidate. The anomaly matrix elements {0\asGG / (4:Tr)\ri') and {0\a gGG / (4:Tr)\G) in Eq. ([50)1 are in reasonable 
agreement with those estimated from the topological susceptibility and quenched lattice calculation. According 
to our analysis, the 77(1405) 77 decay width is 0.6-3 keV, and the leptonic decays 7^(1405) £'^£~ are very 
small. Both predictions can be confronted with future experiments. 
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